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Abstract 

We evaluate the Majorana phases for a general 3 x 3 complex symmetric neutrino mass matrix 
on the basis of Mohapatra-Rodejohann’s phase convention using the three rephasing invariant 
quantities /12, /13 and I2 3 proposed by Sarkar and Singh. We find them interesting as they 
allow us to evaluate each Majorana phase in a model independent way even if one eigenvalue 
is zero. Utilizing the solution of a general complex symmetric mass matrix for eigenvalues and 
mixing angles we determine the Majorana phases for both the hierarchies, normal and inverted, 
taking into account the constraints from neutrino oscillation global fit data as well as bound 
on the sum of the three light neutrino masses (Ejrnj) and the neutrinoless double beta decay 
(PPov) parameter |mn|. This methodology of finding the Majorana phases is applied thereafter 
in some predictive models for both the hierarchical cases (normal and inverted) to evaluate the 
corresponding Majorana phases and it is shown that all the sub cases presented in inverted 
hierarchy section can be realized in a model with texture zeros and scaling ansatz within the 
framework of inverse seesaw although one of the sub case following the normal hierarchy is yet 
to be established. Except the case of quasi degenerate neutrinos, the methodology obtained in 
this work is able to evaluate the corresponding Majorana phases, given any model of neutrino 
masses. 
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1 Introduction 


Apart from hierarchical structure of massive neutrinos a fundamental qualitative nature of these 
elusive particles whether they are Dirac or Majorana type is yet unknown. Neutrinoless double beta 
decay (f3f3o u ) mode [1-12] is able to discriminate between the two different types. Positive evidence 
of the above experimental search will be able to determine the Majorana nature of neutrinos 
assuming the above decay is mediated due to light neutrino. Several f3/3 experiments are ongoing 
and planned. In Ref. [13] a brief discussion about some of the important experiments is presented. 
Among them, EXO-200 [14] experiment puts an upper limit on the relevant neutrino mass matrix 
element |mn| within a range as |mn| < (0.14-0.35 eV). Further, NEXT-100 [15] experiment will 
be able to bring down the above value of the order of 0.1 eV. Thus in an optimistic point of view 
such property of neutrino could be testified by the next generation experiments. However, even if 
it is possible to pin down the value of |mn|, it is still difficult to predict the values of the Majorana 
phases until we can fix the absolute neutrino mass scale. It is shown in Ref. [16] that in addition 
to the /3/3oi, decay experiments, lepton number violating processes in which the Majorana phases 
show up are also corroborative to determine the individual Majorana phases. Another interesting 
physical aspect such as contribution of the Majorana phases to the generation of 6 13 within the 
present 3a range of neutrino oscillation global fit data is also studied in the literature [17]. Ref. [18] 
discusses how to constrain the Majorana phases using the results from cosmology and double beta 
decay. Thus it is worthwhile to study the calculability of the Majorana phases in terms of a general 
neutrino mass matrix (m„) parameters. In the present work we evaluate individual Majorana 
phases in terms of the parameters of a general m u using three rephasing invariants I\ 2 , 113 and I 23 
presented in Ref. [19] on the basis of Mohapatra-Rodejohann’s phase convention [20]. Although 
there are several papers which discusses the general procedure for calculating the Majorana phases, 
motivation behind taking the rephasing invariants is that the methodology we present here is 
capable of calculating the Majorana phase in a model independent way even if one of the eigenvalue 
is zero which is still allowed as far as the present neutrino oscillation global fit data is concerned. 
Moreover as one of the rephasing invariant (I 23 ) is directly proportional to m 3 , therefore it vanishes 
if m 3 = 0 and hence shows a strong dependency of the Majorana phases with the light neutrino 
masses. In the present work we evaluate the Majorana phases for a general complex symmetric 
neutrino mass matrix (mj taking into account the global fit oscillation data and the upper bound 
on the sum of the three light neutrino masses (Ej mi) along with the /3/3o u decay parameter for both 
the hierarchical cases. We then conclude except the case of quasi degeneracy 1 , the methodology 
presented in this work is able to calculate the Majorana phases, given any model of neutrino masses 
and for convenience, we further numerically estimate the ranges of each Majorana phase for both 
types of hierarchies, in the context of a cyclic symmetric model as well as a model with scaling 

1 For the quasi-degenerate case the procedures of calculating the Majorana phases are stated in Sec. 4.3 
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ansatz property. It is also shown that all the sub cases we present in inverted hierarchy section 
of the general discussion can be realized through the choice of a model with scaling ansatz with 
texture zeros within the framework of inverse seesaw while one of the phenomenologically viable 
sub case of the normal hierarchy section is yet to be identified. The plan of the paper is as follows. 

In Section 2 we briefly discuss the basic formalism to set the convention of the Majorana phase 
representation within the framework of neutrino oscillation phenomena. CP violating rephasing 
invariants are presented in Section 3. Section 4 contains explicit calculation of the Majorana phases 
for both types of neutrino mass hierarchies along with phenomenologically viable different sub 
cases. Numerical estimation of the Majorana phases, their connection to the physical observables 
and discussions about their testability for the general case taking into account the constraints from 
the extant data for both types of neutrino mass hierarchies are presented in Section 5. In Section 6 
application of the above methodology in the context of cyclic symmetric and scaling ansatz invariant 
models is presented. Section 7 contains summary of the present work. 


2 Basic formalism 

Experimental observation of neutrino flavour oscillation constitutes a robust evidence in favour 
of nonzero neutrino masses. The flavour transition process is basically a quantum mechanical 
interference phenomena with the explicit relationship between the left handed quantum fields {v a L) 
of the flavour basis and the mass basis {yn) as 


l'aL = ^lUtaiViL (2.1) 

where a(= 1 , 2 ,...., m) corresponds to the flavour index and i{= 1 , 2 ,...., n) implies the mass index 
and the matrix U u is the corresponding neutrino mixing matrix. For three generation of fermions, 
i.e, for n = m = 3, the weak Lagrangian containing charged lepton fields and the neutrino fields 
can be written in the mass basis as 

-c cc =f 2 L L r(u}u:) ai u iL w; + h.c. (2.2) 

where Ui is the unitary mixing matrix in the charged lepton sector. The matrix ujU v is the 
leptonic mixing matrix and is known as the Pontecorvo — Maki — Nakagawa — Sakata mixing 
matrix ( Upmns ) which contains 3 mixing angles and 6 phases in general. It is useful to redefine the 
mixing matrix by absorbing the unphysical phases into the charged lepton fields and the neutrino 
fields (Dirac type). If the neutrinos are Majorana type, they break the global U( 1) symmetry and 
hence, redefinition of the neutrino fields is not possible. Therefore, out of 6 phases 3 unphysical 
phases can be absorbed by redefining only the charged lepton fields and thus the Upmns matrix 


3 


is parametrized as 


Upmns — UckmPm (2-3) 

where Uckm is the usual CKM type matrix and is given by 



^ C12C13 

S12C13 

si 3 e" w \ 


Uckm = 

-S12C23 - Cl 2 S23Sl3e l<5 

C12C23 - Sl2Sl 3 S23e*' 5 

C13S23 

(2.4) 


\ S12S23 - Cl2Sl3C23e l,S 

-C12S23 - Sl2Sl3C23e i<5 

C13C23 / 



where Cij => cos=>■ sin djj and 5 is the Dirac CP phase. Pm is a 3 x 3 diagonal phase matrix 
and following Mohapatra-Rodejohann’s convention [20] it is given by 

P M = (l, e -, e ^+ 5 )) (2.5) 

where a and (3 + 5 are the Majorana phases which do not appear in the neutrino —> neutrino 
oscillation experiments [21,22]. Regarding the structure of Pm matrix we would like to mention 
the following: The advantage of using the above Majorana phase convention is that for m 3 = 0 
it is possible to calculate the single existing Majorana phase a while, for mi = 0, only the phase 
difference a — (/3 + 5) is calculable. The result will be reversed if we utilize the PDG [23] convention. 
Explicitly, with PDG convention, if m 3 = 0, only the phase difference is calculable, however if mi 
is vanishing it is possible to calculate the existing Majorana phase. Based on PDG convention two 
of the authors presented a detailed calculation [24] for both the Majorana phases in context of a 
general rn v . however, if one of the eigenvalue is zero which is still allowed by the present neutrino 
experimental data, it is not possible to calculate individual phases in that case. The above men¬ 
tioned problem is successfully resolved in the present work. 


CP violating effect of Majorana phases in neutrino —> antineutrino oscillation [25-27] and some 
lepton number violating (LNV) processes are studied in detail in Ref. [16]. In this work, using the 
rephasing invariants constructed out of the neutrino mass matrix elements [19] we determine the 
Majorana phases for two different hierarchical cases. 


3 CP violating phase invariants 

Considering neutrinos as the Majorana fermions in extended standard model one can parametrize 
Upmns with the CP violating phases as given in Eqn.(2.3) where we redefine the charged lepton 
fields absorbing the unphysical phases of total mixing matrix U. Hence, in principle the mixing 
matrix U can be defined as 


U — P^Upm NS 


(3.1) 
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where is a 3 x 3 diagonal phase (unphysical) matrix and is given by 


= dja S (e^,e^,e^ 3 ). 


(3.2) 


Now, as the low energy neutrino mass matrix is complex symmetric it can be diagonalized as 


U ] m u U* = d v 


(3.3) 


where 


d v = diag(mi,m2,m^). 


(3.4) 


Substituting Eqn.(3.1) in Eqn.(3.3) we get 


UpMUPMNS 


jt„ 


— di , 


(3.5) 


alternately 

= Up M Nsd v Up MNS . (3-6) 

Thus P^ rotates the mass matrix rri v in phase space. Therefore, the rephasing invariants (remain 
invariant under phase rotation) of m v contain the informations about the CP violating phases. 
It has been shown explicitly in Ref. [19] that for three generations of neutrinos there are three 
independent rephasing invariants and are given by 


I±2 = /?n[m.n?n22Wi2 m 2i] 
h 3 = Im{m22rn^rn* 2Z m% 2 ] 

I 13 = Im[mnm 33 m* 13 m* 31 ] (3.7) 

where m a g is the element of m u at a/3 position with a,/3 = 1,2,3. Now since the invariants of 
Eqn.(3.7) are independent of phase rotation of m u , therefore to evaluate them in terms of mixing 
angles, CP violating phases and the eigenvalues we can rewrite Eqn.(3.6) as 

m u = UpMNsduU pmns (3-8) 


where without any loss of generality we assume (j>i = 0 which corresponds to the structure of as 
Pip = diag{ 1,1,1). Now writing down Eqn.(3.8) explicitly one can find the mass matrix elements 
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as 


mn = cf 2 cf 3 mi + s\ 2 c\ 3 m 2 e 2ia + m 3 s 2 3 e _2 * <5+2 ^ /3+ ‘ 5 ^ (3.9) 

mi2 = Cl 3 {-mi(ci2Sl 2 C23 + c 2 2 si3S 2 3e* <5 ) 

+m 2 e 2ia (ci2Si2C23 - s?2 s i3'S23e i ' S )} + m 3 ci 3 si 3 S23e"'' <5+2 * (/3+<5) (3.10) 

rn 1 3 = Ci 3 {mi(ci 2 Si 2 S23 - C? 2 Sl3C23e Jl5 ) 

-m 2 e 2lQ (ci 2 si2S23 + s?2 s i3C23e i,S )} + m 3 ci 3 si 3 C 2 3 e _ * ,5+2l(/3+<5) (3.11) 

m-n = mi(si 2 c 2 3 + ci 2 S 23 si 3 e * <5 ) 2 

+?n 2 e 2ia (ci2C23 - si2-S23Si3e*' 5 ) 2 + m 3 c 2 3 s| 3 e 2l(/3+,5) (3-12) 

W-23 = mi{ci2Sl2Sl3(c 2 3 - «23) e * 5 + c 12 c 23'S23'Sl3 e2 * <5 — S 12 S 23C23} 

-?n 2 e 2 * a {ci 2 Sl2Sl3(ci 3 - S 2 3 )e tS - sf 2 C 2 3 S 2 3 sj 3 e 215 + C? 2 S 2 3C23} 
+m 3 c 23 S23C? 3 e2l(/3+<5) (3-13) 

m 33 = mi(ci 2 C 2 3 si 3 e t5 - Si 2 s 23 ) 2 

+m 2 e 2ia (si 2 C 23 Si 3 e lS + C 12 S 23) 2 + m 3 c| 3 c 2 3 e 2l(/3+<5) . (3-14) 


It is now straightforward to calculate 1 12 and I 33 using Eqn.(3.9) to Eqn.(3.14). Neglecting terms 
o (4 3) and higher order we obtain /12 and I\ 3 as 

I 12 = Acl 3 lBc23-2s2 3 si 3 {c 2 12 mi<S>i +s\ 2 rri 2 <& 2 }] 

+m|ci2Ci 3 si2C23[—2c 2 2 s 23 cf 3 777.isi 3 ^4i + 2s 23 c 2 3 sf 2 m2Si 3 A 2 ] 

+m 3 Ci2cf 3 Si2C23[si2s| 3 c 13 m l c 12 c 23-43 + 'Sl2 s 23 c 13 m 2Cl2C23 ^4 
-2cf 2 S23"rim2c! 3 Sl3^5 - 2s 2 3"il?Tl2Sl2 c 23 s 13^5 + 2c 2 2 S 2 3S?2 s 13cl 3 ^6 
+2ci 2 s| 3 cl 3 mf si 3 A7 + 4cos(2a)cf 2 s 2 3 ci 3 mi?n2sf 2 si 3 A7 + 2s\ 3 c\ 3 m^s\ 2 si 3 A7\ (3.15) 

Ii3 = Asl 3 [Bs2 3 + 2c 23 si 3 {c 2 l2 mi4>i +s\ 2 m2§2}} 

-mlci2c\ 3 si2S23[-‘2c 2 l2 cl 3 c\ 3 rniSi 3 A 1 + 2cl 3 c 2 l3 m, 2 s 2 l2 si 3 A 2 ] 
+m3Ci2Ci 3 Si2S23[cf 2 c| 3 Ci 3 miSl2S232l 3 + Cl2c| 3 Ci 3 m 2 S? 2 S 2 3^4 
+2Ci 2 C 2 3mim2S 23 Si 3 A 5 + 2c2 3 mim.2s\ 2 sl 3 Si 3 A 3 - 2c 2 2 C23S?2 s 13S 2 3 j4 6 
- 2c i2 c 23 c i3 m i s i3^7 - 4cos(2a)c 2 l2 cl 3 c 2 l3 m l m 2 sl 2 s l3 A 7 - 2c 23 c 2 3 m|si 2 si 3 ^7] (3.16) 

where 


A 

= -ci2Si2mim 2 Ci 3 

(3.17) 

B 

= sin( 2 a)ci 2 Si 2 (m| - m 2 ) 

(3.18) 

$1 

= {sin(2a — 5)mi + sin[<5]m 2 } 

(3.19) 

$2 

= {sin(2a + 5)m 2 — sin [5] mi} 

(3.20) 
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and 


A\ = sin(<5)mi + sin(2a — S)m 2 

A 2 = sin(< 5 )m 2 — sin( 2 a + 8 )m\ 

A 3 = sin 2(/3 + 8 )m\ + 2 sin(2a — 2/3 — 25)m\m2 — sin(4a — 2/3 — 25)m\ 

A 4 = sin 2(a + /3 + 8 )m\ — 2 sin 2(/3 + 5 )m\m 2 — sin(2a — 2/3 — 2(5)m| 

A 5 = sin(2a — 2/3 — 5)m\ + sin(2/3 + S)m^ 

Aq = sin(2/3 + 8 )mf — sin( 2 ai + 2/3 + 5)m\m 2 — sin(4a — 2/3 — 5)mim\ + sin(2a — 2/3 — 5) 

A 7 = sin(2/3 + 5)mi + sin(2o: — 2/3 — 8 )ma- (3.21) 

A careful inspection reveals that the invariants are expressed in a tricky way. To be more precise, 

they are written as 

hj = Ci + 'S13C2 + 11^-3 C3 + ( 3 . 22 ) 

where ‘C % is some parameter dictated by Eqn.(3.15) and Eqn.(3.16). The reason behind such a way 
to write down the invariants are the following: firstly, the popular paradigm in the neutrino mass 
models is to generate vanishing $13 at the leading order and thereafter nonzero value of the same is 
generated by the means of some perturbation to the mass matrix and finally as the oscillation data 
dictates the mass square differences only, there is also a possibility of a vanishing neutrino mass 
(e.g, models with scaling ansatz, Zee-Babu model etc.). Therefore one can see the direct impact of 
their presence or absence in the measures of CP violation. 

Now the remaining invariant (I23) has a special character that it vanishes for m 3 = 0 and it comes 
out as 


I 23 — 777 -3 c 23 s 23 c 13^1 + m 3 c 23 s 23 c 13[2 c 12 c 13 m 2Sl2Sl3(c23 — s 23)-®2 

-2ci2C? 3 miS? 2S i3(4, - s| 3 )-B 3 ] + m 3 C23S23C?3[c? 2 C23"l2 s 23-B4 
+c A V2 C23 rn 1 ml s\ 2 s 23 -65 + c? 2 c\ ? jn\m 2 4 2 s 2 3 B 6 + C23mlsf 2 S23B 7 
-2cf 2 m|si2Si 3 (c2 3 - s\ 3 )B 8 - 4cos(2a)cf 2 mim2S?2' s i3(c| 3 - sl 3 )B 8 ] 

— 2ci2m^sf 2 s 13(c23 — s 23)-®8] (3.23) 
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with 


Bi = sin 2(a- j3- S)cf 2 ci 3 cf 3 m 2 s 2 - sin2(/3 + S)c 23 cf 3 m 1 sf 2 s 23 
= $1 

B 2 = sin(2a — 8)m\ + sin(5)m2 

= - 4>2 

B 3 = sin(5)mi — sin(2a + d)m 2 
B^ = — sin(2a — 2/3 — 26) 

B§ = 2 sin 2(/3 + 6) — sin(4a — 2/3 — 26) 

Bq = — 2sin(2a — /3 — 6) + sin2(2a + j3 + 6) 

Bj = sin2(/3 + <5) 

B% = sin(2/3 + 6)m\ + sin(2a — 2/3 — 6)m 2 . (3.24) 

4 The Majorana phases 

At the outset, first, we would like to mention that the three independent invariants I \ 2 , /13 and 
I 23 stand for the three CP violating phases a, /? + 6 and 6 , however in this section we solve the 
invariants only for the Majorana phases (a, f3 + 6) while the Dirac CP phase 6 is calculable from 
the usual Jarlskog measure of CP violation. Next, for a general m v where all the parameters are 
present and all the eigenvalues and mixing angles are nonzero, all the invariants are independent 
and in principle one can extract the a and f3 + 6 phases without any specific hierarchical assumption 
which is also useful for the quasi-degenerate case. However, the calculation is too cumbersome in 
this general situation. In the present work we consider a simplified approach assuming hierarchical 
structure of neutrino masses and calculate the Majorana phases utilizing the invariants I± 2l /13 and 
I 23 for both, normal and inverted hierarchical cases. 

4.1 Inverted hierarchy (m 2 > m\ >> m 3 ) 

Case I: mi, m 2 , m 3 / 0, $13 / 0: Three independent invariants. 


In this case utilizing Eqn.(3.15) and (3.16)the Majorana phase a comes out as 

/l2S23 + -^13^23 

C 23 S 23 C 13 C 12 S 12 m l m 2 Am| 

where Am|, = m| — m\ and we neglect the terms containing m 3 (m m i n ) in both the invariants (I 12 
and 1 13 ). Another equivalent expression of a can also be obtained from Eqn. (3.9) (neglecting the 
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term containing m 3 s\ 3 ) showing explicit relationship with j3j3 q v decay parameter |mn| as 


1 cos _! f l^nl 2 _ (cf 2 mf + sf 2 ml) \ 

2 \ 2 c^ 2 Si 2 c i 3 r 7 ii m 2 2 c\ 2 s\ 2 m\m 2 J 


(4.2) 


In principle we can use any of the equation (Eqn.(4.1) or Eqn.(4.2)) to find a. The first one depends 
upon the explicit construction of I \2 and /13 in terms of the neutrino mass matrix (m„) elements 
while the second one requires the knowledge of decay parameter |mn|. 

In order to calculate /3 + 5 from Eqn.(3.23) the terms involving si 3 (c 23 — s 23 ) can be neglected. 
Therefore, assuming inverted hierarchy I 23 can be approximated with dominant term as 


h3 = "l2 m -3C23S23C? 3 Ci2C 2 3S23^4 

= -m 2 m 3 c| 3 s 23 c? 3 cf 2 sin(2a - 2 [f3 + 5]). 


Reverting the above equation the Majorana phase (3 + 5 is expressed as 


P 


S = -I sin 


-1 


123 


' n 4 rn 3 c 23 s 23 c 13 c 12 


+ GL. 


(4.3) 


(4.4) 


Case II: mi, m 2 , $i 3 7 ^ 0, m 3 = 0: Two independent invariants. 


In this case utilizing Eqn.(3.15),(3.16) and (3.23) the three rephasing invariants Ii 2 ,Ii 3 and I 23 
come out as 


I 12 = Ac 1 3 [Bc 2 3 - ^S23Sl3{cl 2 mi^l + sl 2 m, 2 ^ 2 }] (4.5) 

= I 12 - 2Acl 3 S23Sl3{c? 2 l«kl + S? 2 m 2 $ 2 } (4.6) 

I\3 = As\ 3 \Bs 2 3 + 2c23Si 3 {c? 2 mi$i + Si 2 m 2 $ 2 }] (4.7) 

= Ji 3 + 2As^ 3 c 23 si3{c? 2 mi<I > i + Si 2 m2 < h2} (4-8) 

I 23 = 0 (4.9) 

where 

I ° 12 = ABc 23 (4.10) 

I \3 = ABs 23 (4-11) 

with A, B already defined in Eqn.(3.17) and (3.18) respectively. As one of the invariant vanishes 
due to the condition m 3 = 0, therefore, the three independent CP phases can not be solved from 
the above invariants and thus the two non zero invariants corresponds to one Majorana phase (a) 
and the Dirac CP phase (4) as ft + 5 vanishes for m 3 = 0. Proceeding as previous we get the 
same expression for the Majorana phase a as given in Eqn.(4.2). Furthermore, solving Eqn.(4.5) 
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to Eqn.(4.7) an equivalent expression of a, same as Eqn.(4.1) is also obtained. 

Case III: mi, m 2 / 0, m 3 , #13 = 0: One independent invariant. 

In this case the invariants given in Eqn.(4.5),(4.7) and (4.9) become 

I vi = I°i2 (4-12) 

h 3 = /i° 3 (4.13) 

and 

h 3 = 0 (4-14) 

It is amply clear that the first two invariants I 12 and I\ 3 are not independent [20] of each other 
and their correlated relationship leads to the estimation of only one Majorana phase a while the 
information about the Dirac CP phase is lost. 

4.2 Normal hierarchy (m 3 >> m 2 > mi) 

Case I: 7771 , 7772 , m 3 / 0, $13 / 0: Three independent invariants. 


In this case since mi = rn mm and 7773 >> 777-2 > mi, we simplify I 12 and I\ 3 as 


I 12 

= 7c sin(2a — 2[/3 + d]) + 77S13S23 sin[d] 


41,3 

= /7sin(2a — 2[/3 + 6]) — 77S13C23 sin[<5] 

(4.15) 

where the parameters k and 77 

are defined through 



2 2 6 3 4 2 

K = -c 12 c 23 c 13 m 2 m 3 .s 12 S 2 3 

(4.16) 


?? = 2ci2C23Ci 3 m2m|sf 2 S23. 

(4.17) 


Now from Eqn.(4.15) we get 

sin( 2 a — 2\/S + - 5 ]) = J 4.1.2 + 4.1.3 1 

= r. (4.18) 

Again due to the hierarchical condition m 3 >> m 2 > mi, I 23 can be approximated as 
I23 - mlc 2 3S23B\ 

= m|c 23 S23[sin(2a - 2 [/3 + 5 ])ci 2 C 2 3 4 3 m 2 S 2 - sin2(/5 + 5)c23c\ 3 m.is\ 2 S23]- (4.19) 
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Inserting Eqn.(4.18) in (4.19) we get 


P + 5 


- sin 
2 


K!2 r _ J 23 

q 9 9 9 

mi "*3 m l c 23 s 23 s 12 c 13 


(4.20) 


where cti 2 => cot 0 i 2 - 

It is now straight forward to calculate the other Majorana phase a from Eqn.(4.18) and it comes 
out as 


sin ” 1 T + 2(/5 + 5) 


(4.21) 


Case II: m 2 , m 3 , 6*13 7 ^ 0, mi = 0: Two independent invariants. 


In this case neglecting terms like sf 3 and si 3 (ci 3 — s 23 ) in I 23 only the Majorana phase difference 
(a — [p + d]) is calculable and is given by 


a-[P + 5} 


- sin 
2 


h^mls \ 2 \ 
— nm 3 J 


(4.22) 


along with an explicit relationship between the three invariants as 


I 23 


m% 


C23I12 + S 23 Il 3 


m 2 sin 4 0 i 2 


(4.23) 


Therefore, essentially we get two independent invariants corresponding to the Majorana phase dif¬ 
ference and the Dirac CP phase. 


Case III: "* 2 , 777.3 / 0, mi, 0 i 3 = 0: One independent invariant. 

In such a condition the three invariants are coming out in a correlated manner as 

-612 = 

I 23 = 

and in this case only independent invariant I\ 2 is connected to the Majorana phase difference 
(a — [p + 5]). 


«sin2(a — [P + <5]) 

- sin 2 (a - [P + 5})cj 2 4 3 sl 3 cl 3 s\ 2 mlm 3 

Il3 

sin 2 (a - [P + 5})cl 2 cl 3 sl 3 ci 3 m 2 inl 


mi 

2*12 1 

| -6l2 

m 2 ) 


2 \ 


m 3 4 1 

2 *12 1 

| ^13 

m 2 ) 



(4.24) 
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4.3 Quasi-degenerate case 


Although in the present work we are not discussing the quasi-degenerate case which is relevant in the 
cosmological context [28], however, one can calculate the Majorana phases in a model independent 
way by directly solving the invariants as mentioned at the beginning of Sec. 4. To be precise, 
using Eqn. (3.15), (3.16) and Eqn. (3.23) one can extract all the CP violating phases without 
any hierarchical assumption. However, the calculation is tedious and will be studied elsewhere. 
Another alternative way is to follow the calculations presented in Ref. [24]. However, in that case 
the phase convention is different. Utilizing the phase convention presented in this work one can 
calculate all the phases in the second method also. 


5 Numerical estimation 

5.1 Parametrization, diagonalization and the ranges of the Majorana phases 

A general solution for a three generation complex symmetric Majorana mass matrix is given in 
Ref. [24]. In order to estimate the Majorana phases obtained in the present work we utilize the 
expressions of the three eigenvalues and the three mixing angles. We also use the global fit data of 
neutrino oscillation experiments shown in Table 1 and the upper limits on the sum of the neutrino 
masses (E,mi(= mi + m 2 + m 3 ) < 0.23 eV) [29] and the /3/?o^ parameter (|mn| < 0.35 eV) [14] to 
obtain model independent ranges of the Majorana phases. 


Table 1: Input experimental values [30] 


Quantity 

3d ranges 

|Am|i| (N) 

2.31< Am3 1 (10 3 eU -2 ) < 2.74 

|Am|i| (I) 

2 . 21 < Am| 1 (10 3 eU" 2 ) < 2.64 

Am], 

7.21< Am2 1 (10 5 eU -2 ) < 8.20 

012 

31.3° < 0\2 < 37.46° 

023 

36.86° < 023 < 55.55° 

013 

7.49° < 6 >i 3 < 10.46° 


We consider a most general 3x3 complex symmetric neutrino mass matrix m u as 



(P 

Q 

R 

m v = 

Q 

s 

T 


\R 

T 

V 
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(5.1) 











with all parameter complex and can be parametrized as 



( 1 

xe lolx 

ye loly ^ 



m v = moe*“ m 

xe iax 

ze ia * 

we lolw 


(5.2) 


\ye lcXy 

we lolw 

ve lolv ) 



with the definition of the parameters 






P = m 0 e ia ™,Q/P = xe iotx , R/P = ye™* 

n 

--ze ia % 

T/P = we ia ™,V/P = ve iav . 

(5.3) 


We can now give a phase rotation to the matrix of Eqn.(5.2) by a diagonal phase matrix K = 
diagte 1 ^ 1 , e*^ 2 , eifa) as 


f ' I' 

m„ = K m u K 


and consequently the rotated matrix comes out with 9 parameters as 


/I 


m u = m 0 


X 


y 




x ze llL1 we 
\y we lii2 ve 1 ^ 3 


where x,y,z,w are the real parameters and Hi, ^ 2 , ^ 3 , <t> 1 , 4> 2 , 4>3 are defined as 


(5.4) 


(5.5) 


11 1 — ot z ‘2otxi ^2 — ®-x ^3 — (5.6) 

and 

*1 = = -{a x - ^),h = -K - ^)- (5-7) 

Now using Eqn.(3.7) we can explicitly calculate the rephasing invariants in terms of the elements 
of m' u . It is to be noted, that in the general case the number of parameters are 9 and we have 
only 7 experimental inputs. However, among the 9 parameters there are three angle parameters 
(Hi, H 2 and H 3 ). We set the values of these angle parameters in an arbitrary manner within the 
range 0 — 27r and vary the other parameters in a wide range to estimate the overall ranges of the 
Majorana phases which are depicted in figure 1. We first constrain the rephasing invariants which 
in turn generate the correlated plot of the Majorana phases. The correlation between the phases 
are the consequences of Eqn.(4.21) and Eqn.(4.4) respectively. 

Upon numerical estimation, the model independent ranges for a and f3 + 5 come out as —90° < 
a < 90° and —71° < /3 + S < 71° for normal hierarchy (mi / 0, 6*13 7 ^ 0) and —45° < a < 45°, 
—70° < (3 + 5 < 70° for inverted hierarchy (m 3 7 ^ 0, $13 7 ^ 0) and are shown explicitly in figure 1. 
For m 3 = 0 case, the range of a is obtained as —45° < a < 45° and for case mi = 0, the phase dif¬ 
ference is constrained as —82° < a— [/3+<5] < 82°. We also present the parameter ranges in Table 2. 
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Figure 1: Plots of the Majorana phases (a vs /? + <5) for normal (left) and inverted (right) hierarchies . 


Table 2: Parameter ranges for a phenomenologically viable m u . 


Hierarchies j. 

m 0 

X 

y 

z 

w 

V 

Normal 

A 0 

2.4 x 10" 4 - 
1.7 x 1(T 3 

0.15 < x < 

3.7 

0.15 < y < 

4.6 

0.14 < z < 9.5 

0.1 < w < 

8.6 

0.13 < v < 

8.4 

Normal 

m\ = 0 

2x 10 -4 —1.2x 
1CT 3 

0.1 < x < 3.2 

0.14 < y < 

4.7 

0.09 < z < 7.5 

0.09 < w < 

8 

0.11 < v < 

8.1 

Inverted 

m 3 a 0 

1.2 x 10" 4 - 
1.8 x 1(T 3 

0.5 < x < 3.5 

0.5 < y < 

3.47 

0.1 < z < 2.6 

0 < w < 1.8 

0 < v < 2.4 

Inverted 

m 3 = 0 

1.1 x 10“ 4 - 
1.4 x 1(T 3 

0.1 < x < 3 

0.2 < y < 3.4 

0 < z < 2.4 

0 < w < 1.7 

0 < v < 2.4 


5.2 Connection to the physical observables and future of the Majorana phases 

As previously mentioned, unlike the Dirac CP phase 5, the Majorana phases do not appear in the 
neutrino—^neutrino oscillation. Therefore, a natural question arises how and where these phases can 
be measured. As a direct detection, in Ref. [25] Xing suggested a thought experiment (neutrino —> 
antineutrino oscillation) in which he pointed out these phases may appear in the probability expres¬ 
sion of the flavour oscillation and thus also in the expression of the CP asymmetry parameter A a p 
which is the measure of CP violation. However, this kind of experiment is purely academic at this 
moment and practically difficult to design as the oscillation probability is highly suppressed by the 
factor m 2 /E 2 , where rn t is the mass of the light neutrino and E is the beam energy. Now consider¬ 
ing E ~ MeV and the masses of the neutrinos to be less than 1 eV, one can calculate mf/E 2 to be 
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e>(10 _12 ). To improve rrii/E, a novel suggestion [16,25] is to lower the value of E, however, in that 
case the estimated size of the base line length and the detector are beyond the reach of the present 
experimental facilities. But as an optimistic point of view we expect these kind of experiments 
will be designed in future and thus the prediction of the Majorana phases will be tested. Beside 
neutrino —> antineutrino oscillation there are several LNV processes like (5(3o„ decay, A ++ —» l+lp 
(in Type II seesaw model) [16] etc., which play a crucial role for the indirect measurement of the 
Majorana phases. Now coming into our work, we present a table in the appendix which shows 
the ranges of the obtained Majorana phases for some typical values of |mn| and for convenience, 
in figure 2 we present variation of the Majorana phases with |mn| for the best fit value of Ara^ 
and taking all the other constraints in their 3a ranges for both the hierarchies. We would like to 
mention that even if we take the 3a range of Am|i, over all ranges of the Majorana phases do not 
differ much, however, unlike the plots of figure 2, the plots in that case become more wider for the 
higher values of |mn| (> 0.08 eV). Although, the present experimental upper bound on |mn| is 
0.35 eV, NEXT will be able to bring down the value to 0.1 eV and thus the approximate ranges of 
the Majorana phases can be predicted. 

Thus far we have estimated the Majorana phases in a general context. Latter, we apply the 
expressions obtained for a and (5 + 5 for few testable flavour models (models with lesser number of 
parameters) as an application of the general result although our analysis is true for any hierarchical 
model of neutrino masses. 


6 Some testable flavour models 

The reason we discuss this section is to make certain whether the results obtain in the general case 
are consistent with the other models or not. Moreover the models with certain flavour symmetries 
are highly predictive in nature. Therefore, precise measurement of the CP violating phases may 
act as the important tools to verify the testability of the flavour models [31]. In inverted hierarchy 
section we present a model with scaling ansatz and texture zeros within the framework of inverse 
seesaw through which all the sub cases presented in Sec. 4.1 can be realized while in the normal 
hierarchy section we present a model with cyclic symmetry within the framework of Type I seesaw. 
Obviously the choices are for illustration. One can also consider inverse or linear seesaw for normal 
hierarchy [32,33] and Type I seesaw for inverted hierarchy [34], In principle one can use the 
technique in any hierarchical flavour models. The above numerical results are obtained for the 
general m„ where all the 9 independent parameters are present. However, as previously said, one 
can reduce the number of parameters by invoking some symmetry or ansatz in the Lagrangian 
which is more predictive in nature and thus testable in the experiments. In this section we provide 
applications of the general results in few typical cases for both the hierarchies, normal and inverted. 
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Figure 2: Plots of the Majorana phases (a, /3 + 5) vs |toh| for normal (left) and inverted (right) hierarchies 
for best fit values of Am^. 

6.1 Normal hierarchy 

In this case we explore a model that corresponds to Case I of the normal hierarchical scenario 
mentioned in section (4.2). The model is based on cyclic symmetry with Type I seesaw mechanism 
to accommodate the neutrino oscillation data. In the fundamental level the symmetry exists in the 
neutrino sector of the Lagrangian and due to the symmetry a degeneracy in masses occurs removal 
of which therefore requires breaking of the symmetry. It is shown that a minimal breaking in the 
Majorana mass matrix is sufficient to explain the extant data. In this model the low energy broken 
symmetric mass matrix m v {= —moM^rrij^) originated from Type I seesaw mechanism is given 
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fyi 

2/2 2/3 \ 


m D = 

2/3 

yi 2/2 1 , 
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\V2 

2/3 yij 


= diag(m + £ 1 , m + £ 2 , m) 

(6.3) 


with ei and £2 as the breaking parameters and 


mo 


= -^,pe iQ = ae ip = 2/2 . 


qe H = 




(6.4) 
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Figure 3: Correlation of a vs /3 + <5. 


For numerical analysis we choose the mass scale of Mr to be of the order of 10 15 GeV and rno 
to be at electroweak scale. Further redefining the breaking parameters as e[ = — and e' 2 = ^ we 
allow them to vary as —0.1 < £^,£2 < 0.1 to keep the breaking effect small. We then constrain 
the parameter spaces taking into account the 3c ranges of neutrino oscillation global fit data and 
explicitly evaluate both the Majorana phases. From figure 3 the ranges read as —77.2° < a < 76.7° 
and —45.3° < j3 + 5 < 45.5°. Note that the ranges of both the phases are embedded within the 
values obtained for the general case. Similar to the general case, in figure 4 we also present the 


17 
























variation of the Majorana phases with the parameter. One can see the upper limit of |mn| is 
~ 0.07 eV which is well within the reach of the future planned experiments. 

As the model consists of lesser number of parameters, we also expect a significant correlation 



0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 


|mii | (eV) |mn|(eV) 


Figure 4: Variation of a and j3 + S with |toh| for cyclic symmetric case (normal hierarchy). 


between the phase invariants and are depicted in figure 5. 




1 12 (in (eV) 4 ) 1 13 (in (eV) 4 ) 

Figure 5: Correlated plots of the rephasing invariants, I\2 vs I 23 (left) and /13 vs I23 (right). 
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6.2 Inverted hierarchy 


In this case , we explore a model based on scaling ansatz with inverse seesaw mechanism [32,35-46] 
. In this mechanism m u is given by 

m v = 

(6.5) 


where mr> is the usual Dirac type matrix and the other two matrices /x (Majorana type) and Mrs 
(Dirac type) arise due to the interaction between the additional singlet fermion and right handed 
neutrino considered in this type of seesaw mechanism. To further reduce the number of parameters 
texture zeros [33,34,47-74] are assumed in the constituent uir and /x matrices. Scaling ansatz 
invariance dictates m 3 = 0 and @13 = 0 and this case corresponds to Case III of Sec. 4.1. Thus to 
generate non zero #13 breaking of the ansatz is necessary. Incorporating breaking in mo through a 
small parameter e, there are two different phenomenologically survived textures which are given by 
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where all the parameters are complex [74]. In both the cases $13 7 ^ 0 however, m 3 = 0 due to 
singular nature of p matrix and this case corresponds to Case II of Sec. 4.1. 

We further consider the most general version of the above case through the breaking of the ansatz 
in both mr> and p matrices through two small parameters e and e' respectively and the neutrino 
mass matrix ml comes out as 


ml 


/ 1 

k\p 

P \ 


(0 

0 

° 

k\p 

k 2 (q 2 e l9 + p 2 ) 

ki(q 2 e zd + p 2 ) 

+ m 0 e 

0 

2 k 2 q 2 e l9 

k\q 2 e l9 

V P 

ki(q 2 e ld + p 2 ) 

(q 2 e l 9 +p 2 ) ) 



to 

Ct) 

Qs 

0 ) 


+moe 


/ 0 kip p 1 
k±p 0 0 

\ p 0 0 ; 


( 6 . 8 ) 


19 









and in this situation both 613 and m 3 are nonzero corresponding to Case I of Sec. 4.1. Thus the 
whole inverted hierarchical sector is generated through the choice of the above model. 

Now, with the explicit construction of rephasing invariants we calculate the Majorana phases in 
each case. Interestingly, for all the cases, the value of Jcp comes out very small due to smallness of 
the Dirac CP phase 5, or more precisely, due to almost real nature of the mass matrices. Therefore, 
such typical nature of the mass matrices also constrains the Majorana phases approximately as 
— 1.2° < a < 0.8° for the first two matrices (ml and ml) and —0.17° < a < 0.17°,—1.5° < 
P + 5 < 1.5° for the matrix ml along with an approximate range of PPo u parameter |mn| as 0.01 
eV < |mn| < 0.0148 eV and 0.01 eV < |mn| < 0.0152 eV respectively. For illustration, in figure 
6 we plot a and [3 + 8 with |mn| for ml- For other two matrices (ml and ml) the variations of 
a with |mn| are almost same as that of the extreme left plot of the lower panel of figure 6 . The 
model is highly predictive and hence, if significant CP violation is observed, the model will be ruled 
out. We plot the correlation between the invariants in figure 7. 

Finally, we summarize our results in Table 3 that shows the ranges of the Majorana phases for all 
the cases. 

Table 3: Summary of the numerical results. 


Hierarchies —> 

Particular cases —> 

General case 

Cyclic symmetry 

Scaling ansatz 

Normal 

Inverted 

Normal 

Inverted 

Inverted 

m\ / 0 

#13 / 0 

m 3 yl 0 

013 / 0 

mi ^ 0 0i 3 / 0 

m 3 / 0 

013 7 - 0 

m 3 = 0 

013 ^ 0 

a ( deg.) 

-90 - 90 

-45 - 45 

-77.2 - 76.7 

-0.17-0.17 

-1.2 -0.8 

P + 8 ( deg.) 

-71 - 71 

-70 - 70 

-45.3- 45.5 

-1.5 - 1.5 

absent 


7 Summary 

In the present work we calculate the Majorana phases of a general complex symmetric 3x3 
neutrino mass matrix utilizing the three rephasing invariant quantities J 12 , I 13 and I 23 proposed by 
Sarkar and Singh for both the hierarchical structures of neutrinos using Mohapatra-Rodejohann’s 
phase convention. Motivation behind the usage of the invariants to calculate the Majorana phases 
is that such methodology enables us to evaluate the existing Majorana phase even if one of the 
eigenvalue (m 3 ) is zero in a model independent way. However, if m± = 0, this methodology will 
only enable us to calculate the difference of the Majorana phases therefore it is needed to change 
the phase convention in that case. After the presentation of the generalized prescription, we further 
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Figure 6: Correlation plot of a vs j3 + 6 (upper panel) and variation of a and /3 + 6 with |mn| for inverted 
hierarchy : scaling ansatz case (lower panel). 


present the maximal allowed ranges of the Majorana phases in general context for both the viable 
hierarchical structures of neutrino masses and address our methodology to be true for any model 
except the case of quasi degeneracy. We then talk about the connection of the Majorana phases with 
physical observables like parameter \mu\ and the branching ratios of charged Higgs (A ++ ) 
decay where the phases show up. As a direct measurement of the Majorana phases we give the 
example of neutrino —> antineutrino oscillation which is a thought experiment right now, however, 
well studied in literature. Although the presented methodology can be used in any hierarchical 
models in neutrino physics, after discussing the general case we further exemplify our methodology 
in few typical testable models (models with lesser number of parameters) leading to normal and 
inverted hierarchy and with their significant predictions on |mn| and the Majorana phases. For 
normal hierarchical case we give an example of a model based on cyclic symmetry with Type I 
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Figure 7: Correlated plots of the rephasing invariants, I \2 vs J 23 (left) and /13 vs I 23 (right). 


seesaw mechanism. We estimate the Majorana phases for the broken symmetric case, since cyclic 
symmetry dictates a degeneracy in the mass eigenvalues. As an example of inverted hierarchy, we 
cite a model with high predictability and comprised of scaling ansatz, texture zeros and inverse 
seesaw mechanism. It is seen that all the sub cases of inverted hierarchy mentioned in Section 4.1 
can be obtained depending upon the scheme of incorporation of ansatz breaking mechanism while 
a phenomenologically viable sub case (mi = 0) of the normal hierarchy is yet to be established 
through the choice of a suitable model. 
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A Appendix 


Table 4: Majorana phases for the general rn v for some typical values of |mn|. 


Normal Hierarchy 

\m xl \ (eV)-» 

0.001 

0.005 

0.01 

0.05 

0.1 

0.12 

0.15 

0.2 

0.25 

0.30 

a (deg.) 

-90- 

90 

-90- 

90 

-87- 

87 

-82- 

82 

-80- 

80 

-77- 

77 

-76- 

76 

-75- 

75 

-75- 

75 

-74- 

74 

P + 5 

(deg.) 

-71- 

71 

-71- 

71 

-71- 

71 

-70- 

70 

-70- 

70 

-69- 

69 

-67- 

67 

-es- 

es 

-65- 

65 

-64- 

64 

Inverted Hierarchy 

|mn| (eV)—>• 

0.01 

0.03 

0.05 

0.08 

0.10 

0.12 

0.15 

0.2 

0.25 

0.30 

a (deg.) 

-70- 

70 

-70- 

70 

-70- 

70 

-69- 

69 

-69- 

69 

-69- 

69 

-es- 

es 

-65- 

65 

-es- 

es 

-63- 

63 

P + 5 

(deg.) 

-46- 

46 

-46- 

46 

-46- 

46 

-45- 

45 

-45- 

45 

-44- 

44 

-43- 

43 

-42- 

42 

-42- 

42 

-40- 

40 
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